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We use anhydrous ferric chloride (FeCk) to intercalate graphite flakes consisting of 2 to 4 graphene 
layers and to dope graphene monolayers. The intercalant, staging, stability and doping of the 
resulting intercalation compounds (ICs) are characterized by Raman scattering. The G peak of 
monolayer graphene heavily-doped by FeCl3 upshifts to^ 1627cm - 1 . 2-4 layer ICs have similar 
upshifts, and a Lorentzian lineshape for the 2D band, indicating that each layer behaves as a 
decoupled heavily doped monolayer. By performing Raman measurement at different excitation 
energies we show that, for a given doping level, the 2D peak can be suppressed by Pauli blocking 
for laser energy below the doping level. Thus, multi-wavelength Raman spectroscopy allows a direct 
evaluation of the Fermi level, complementary to that derived by Raman measurements at excitation 
energies higher than the doping level. We estimate a Fermi level shift of~0.9eV. These ICs are ideal 
test-beds for the physical and chemical properties of heavily-doped graphenes. 



Graphite intercalation compounds (GICs) are host- 
guest systems [l|, Q. They have been intensively stud- 
ied because of novel features in their structures, elec- 
tronic and optical properties [1-7]. They are promis- 
ing for electrodes, conductors, superconductors, cata- 
lysts, hydrogen storage materials, batteries, displays, 
polarizers pHa, Q. Since the first synthesis of GICs in 
1841@, hundreds of GICs have been produced with a 
variety of donor and acceptors intercalantsp], Q. GICs 
featuring unique stackingsequences (staging) have also 
been intensively studied [l|-l4| . In stage 1 GICs, each 
graphene layer is sandwiched by two intercalant layers [l|- 
[j, However, it is difficult to manipulate and pro- 
cess traditional GICs into nanoelectronic devices due 
to their thickness [H, Q. Graphene ha s g reat potential 
in nanoelectronics and optoelectronics [Tol [Til ] . By in- 
tercalating graphite flakes just a few layers thick, one 
can combine the physical and chemical properties of 
GICs with those of few-layer graphenes (FLG) and open 
new opportunities for applications in nanoelectronics [1 2\- 
2l|. There is also great interest on the transport 
properties of graphene at high carrier density, both 
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[19[. By means of an electrolytic gate, Refs llil 1311 18l 
doped graphene up to~4.5x 10 13 cm -2 . Refll9l used 
an ionic-liquid gate to achieve a carrier density higher 
than 10 14 cm -2 . Refll7l achieved 4xl0 14 cm -2 for hole 
and electron doping by means of solid polymer elec- 
trolyte. We note that in donor- type graphite interca- 
lation compound, such as KCg, the electron density can 
reach up to~4.8x 10 14 cm -2 , corresponding to Fermi shift 
There is thus scope for using a similar 
approach to achieve graphene doped at levels higher than 
those reported in RefsJa, [l2, 13.Tl7-19[ 




Here, we use FeCl3 to intercalate FLG and dope sin- 
gle layer graphene (SLG). Raman spectroscopy at several 
wavelengths shows the formation of acceptor-type Stage- 



FIG. 1. a) Optical image of a pristine 1-4L flake, b) AFM 
image of the region indicated by a solid square in (a), c) 
Optical image after FeCl3 doping/intercalation. The scale 
bar in (a,c) is 4/xm. 



1 GICs. We estimate a Fermi shift of~0.9eV, correspond- 
ing to a fractional charge transfer of~l/6.6 = 0.152 holes 
per carbon atom, i.e., a hole density of~5.8x 10 14 cm -2 . 

Graphite flakes consisting of 1 to 4 layers (L) are ob- 
tained by micro-mechanical cleavage of natural graphite 
on a Si+300nm Si02 substrate 1221. The number of layers 
is identified by optical contrast [23l Q and atomic force 
microscopy (AFM) [22^ Figs. la, b. The Raman spec- 
tra are taken at room temperature using a Jobin-Yvon 
HR800 system with a~1.2cm _1 spectral resolution. In- 
tercalation is performed following the vapor transport 
method commonly used in GICs, as discussed, e.g., in 
Refill. The intended intercalant and the flakes are po- 
sitioned in different zones in a glass ampoule, as shown 
in Fig. 5b of Refjl|, then pumped to 1.5x 10 _4 Torr. The 
reaction is carried out at 613K for six and thirty hours, 
for FLG and bulk graphite, respectively. A longer time is 
needed to reach stage 1 GIC for bulk graphite, due to the 
sample size, both in spatial extent and thickness. Fig.lc 
is an optical micrograph of a representative intercalated 
flake. The number of layers can still be identified, and 
their optical contrast is higher than prior to intercalation. 

In GICs only a few layers thick it is difficult to apply 



2 



X-ray diffraction techniques, unlike for bulk GICs staging 
determination (H-Q, because of the small flake thickness 
and the resulting substrate effects. Raman scattering was 
used to distinguish the intercalation and adsorption be- 
havior of Bromine (Br2), Iodine (12), FeCl3 and sulfuric 
acid (25l-[27| . In principle, for stage-1 GICs, a single G 
peak is expected [iMfi- However, multiple G peaks were 
reported in recent works [26). 

The Raman spectrum of graphene consists of a set of 
distinct peaks. The G and D appear around 1580 and 
1350cm -1 , respectively. The G peak corresponds to the 
Fj2 9 phonon at the Brillouin zone center. The D peak is 
due to the breathing modes of six-atom rings and requires 
a defect for its activation [28- 301. It comes from TO 
phonons around the K point [HI, I23, [Hj , is active by dou- 
ble resonance (DR) [30] , and strongly dispersive with ex- 
citation energy due to a Kohn Anomaly at K[32]. DR can 
also happen intra-valley, i.e. connecting two points be- 
longing to the same cone around K (or K^jH, This 
gives the so-called D' peak, which is at^ 1620cm -1 in de- 
fected graphite measured at 514nm[35]. The 2D peak is 
the second order of the D peak. This is a single peak in 
SLG, whereas it splits in four in bilayer graphene (BLG), 
reflecting the evolution of the band structure [36|. Ra- 
man spectroscopy allows monitoring of doping, defects, 
strain, disorder, chemical modifications, edges, and rela- 
tive orientation of the graphene layers [12l ll 
I45I ]. Each Raman peak is characterized by its position, 
width, height, and area. The frequency-integrated area 
under each peak represents the probability of the whole 
process [46|. We thus consider both area, A(2D)/A(G), 
and height, I(2D)/I(G), ratios. 

The G peak position, Pos(G), has been widely used 
to identify s t aging [l|-[4|. In graphene, the shift of the 
Fermi energy has two major effects: (1) a change of the 
equilibrium lattice parameter with a consequent stiffen- 
ing/softening of the phonons [T2I. [l3l. l47| . and (2) the on- 
set of effects beyond the adiabatic Born-Oppenheimer 
approximation, that mo dify the phonon dispersion close 
to the Kohn anomalies 0, [H SI, El S, (H, Hsj. (2) 
always results in a G upshift, both for electron and hole 
doping 32l. 137^ while (1) gives an upshift for p doping 
and a downshift for n doping. Thus for low doping 
levels, below^3 x 10 13 cm~ 2 , both n and p doping result 
" EE S, EI 5ijj|, but for levels 
the n doped case reverses [49]. In 



in G peak upshifts 12, 
above~3xl0 13 cm -2 

fact, for n doping~10 14 cm -2 the G shift would revert to 
zero [49| . Doping also significantly affects the 2D peak in- 
tensity. Ref. |46| reported that, for excitation energy well 
above the Fermi energy, this is linked to A(2D)/A(G) as: 



1 A(G) _ 0.26 

A(2D) ~ Je-ph 
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FIG. 2. Raman spectra of pristine (dashed lines) and 
doped/intercalated (solid lines) 1-4L flakes, measured for 
532nm excitation, a) Low frequency region, (b) D- and G- 
region. c) 2D-region. In (a) we also report the Raman spec- 
trum of bulk FeCl3 (dotted gray line) for comparison. Vertical 
dotted lines are guides to the eye 



28|, |32|, |36|- charge, electron velocity and dielectric constant 46]. 



FeCl3 has eight Raman- active modes 4A g +4E g 
(2Ai p +2A 2 p+4E p )[3]. Six Raman modes (3A^+3E^) 
were thus far experimentally measured [3]. When FeCl3 
was used as intercalant in stage 1 GICs, only four Ra- 
man modes (2Ai^+2E^) were observed [3]. The other 
two A^ and E^ modes at~164 and 354cm -1 are prob- 
ably too weak to be observed in GICs. Indeed, even in 
crystal FeCl3, those modes are very weak[3j. Fig. 2a 
shows that, after doping by FeCl3, three Raman modes 
from FeCl3 are observed in the low frequency region: 
Ai3~136cm -1 , E -179cm -1 and Ai -287cm -1 . We 



where ^e-ph, Ep, e, vp and e are the scattering rate 
due to the emission of phonons, Fermi energy, electron 



cannot detect the other E^ mode~93cm -1 , since this 
is too weak to be distinguished from the background. 
These peaks positions are very close to those previously 
observed for FeCl3 intercalated Stage-1 GICs[l|-[3[, but 
differ from bulk crystalline FeCl3, whose spectrum is also 
shown in Fig. 2a for comparison: the A\ g (at 136cm -1 ) 
and Eg mode (at 179cm -1 ) have a~4cm -1 downshift, 
while the A\ g (at 287cm -1 ) upshifts~3cm -1 . In bulk 
FeCl3, the iron layer is sandwiched by two chlorine lay- 
ers, as shown in Fig.l of Ref|3|. When intercalation hap- 
pens, CI atoms simultaneously occupy preferred sites as- 
sociated with the graphene lattice, which results in the 
loss of the CI atoms long-range two dimensional order, 
because their in-plane structure is incommensurate with 
the graphene host lattice [l], Fe atoms however retain 
the long-range order as in crystal FeCl 3 [E H- This re- 
sults in a~3° relative rotation of the CI layers above and 
below the Fe layer, and a difference of the Raman modes 
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FIG. 3. a) G and b) 2D band of Stage- 1 flakes with 2/3L, and 
doped SLG, measured at 532nm for samples kept in the sealed 
quartz tube used for intercalation/doping, c) Schematic dia- 
grams of FeCl3 doped/intercalated 1-3L flakes [1] 



of intercalated FeCl3 compared to bulkFeCls [l|, Q . FeCl3 
modes are not observed in the FeCl3-doped SLG due to 
the very low density of FeCl3, compared to fully- FeCl3- 
intercalated GICs, due to de-adsorption. 

Figure 2 plots the Raman spectra of 1-4L flakes mea- 
sured at 532nm before (dashed lines) and after (solid 
lines) FeCl3 intercalation. The pristine samples have 
the characteristic features of mono- and multi-layer 



graphene 36, 39]. The Pos(G) shift of doped/intercalated 
samples compared to pristine ones in Fig. 2b is a signa- 
ture of doping. The blueshift is smaller in SLG com- 
pared to FLG. This indicates fewer FeCl3 molecules on 
SLG relative to 2-4L flakes. The SLG 2D band in Fig.2c 
upshifts~28cm _1 , typical of hole- doping (l2i [l3| . while 
I(2D)/I(G) and A(2D)/A(G) decrease-61% and 53% rel- 
ative to those prior to doping. From RefsHH and H, we 
estimate the Fermi shift for SLG to be~0.4eV. 

The 2D lineshape for the 2-4L flakes after FeCl3 inter- 
calation changes significantly, as shown in Fig. 2c, from 
multiple peaks to a single Lorentzian. This is an indica- 
tion of electronic decoupling of the layers (3fj l39| . Note 
that the residual presence of any non-intercalated BLG 
or FLG would give residual multiple 2D bands. In pris- 
tine Bernal-stacked graphite, the interlayer distance is 
When FeCl3 molecules are intercalated, the 
distance increases to 9.37iL[l|-l4(. As a result, the inter- 
layer interaction significantly decreases [H 0, |fj. There- 
fore, the single Lorentzian 2D peak indicates SLG be- 
tween two intercalant layers, each hole-doped. 

As shown in Fig. 2b, the G bands of SLG and BLG 
are~1605 and~1615cm _1 , moving to~1625cm _1 for 3L 
and 4L. The latter value is close to that previously re- 
ported in FeCl 3 Stage-1 GICs (-1626cm- 1 ) (JQ. For 
3L and 4L, an additional sideband appears~1618cm _1 , 
with almost equal width to the main peak. One might 
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FIG. 4. Raman spectra for 532nm excitation energy of as- 
prepared 2L-GICs and after 12 hours, and 1, 4, 6, 8 months: 
(a) Low frequency region; (b) D and G band region; (c) 2D 
band region. A pristine BLG is also included as reference. 
Vertical dotted lines are guide to the eyes. 



assign it as the D' peak activated by defects. However, 
this is not the case for our samples since the D peak 
is unobservable for all the layers, before and after dop- 
ing/intercalation, Fig. 2b. This band is thus another G 
peak, resulting from non-uniform intercalation, due to 
de-adsorptionl of FeCl3 when the flakes are exposed to 
air. This would imply that the shoulder is due to the 
upper or lower layers. The thicker the flakes, the less the 
top and bottom layers will contribute to the overall inten- 
sity of the measured Raman spectrum. Indeed, 4L flakes 
have weaker relative intensity of this shoulder compared 
to the main G peak, than 3L flakes. No such shoulder is 
expected in SLG, where the lower Pos(G) is an indication 
of lower coverage. We note that no shoulder is observed 
for the doped/intercalated BLG, indicative of homoge- 
nous doping. Let us consider the upper and lower layers 
of intercalated 2-4L flakes. If FeCl3 is only present on 
one side of these layers, we expect the amount of charge 
transfer to reach at most that of Stage-2 GICs. In this 
case Pos(G) can only shift to~1612cm _1 , i.e. Pos(G) 
of FeCl 3 -intercalated Stage-2 GICs[l|-[3. Pos(G) of our 
FeCl3-doped SLG is very close to that reported in Reflirj 
where double G peaks, at~1612 and~1623cm _1 were also 
observed for 3 and 4L samples. Refl26l explained this by 
arguing that FeCl3 did not adsorb on the top and bottom 
of their flakes. We note that Pos(G) of our intercalated 
BLG after air exposure (~1615cm _1 ) and of the lower en- 
ergy G band in 3-4L FeCl3-intercalated flakes are higher 
than previously reported for FeCl3 intercalated Stage 2 
GICs[l|-[4(. This means that the upper and lower sur- 
faces of our FeCl3-intercalated 2-4L are doped, but not 
with full coverage, due to FeCl3 de-adsorption, unlike the 
interior layers. The lower Pos(G) in FeCl3-intercalated 
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FIG. 5. I(2D)/I(G) and A(2D)/A(G) as a function of Pos(G). 
The blue and red dashed lines are guide to the eyes 



BLG compared to 3 and 4L is expected, since a BLG 
has no fully enclosed SLG, thus de- adsorption effects are 
expected to be stronger. 

If the samples are quickly removed from the sealed 
tube and cooled to room temperature in air, their Raman 
spectra are very similar to those previously reported for 
Stage- 1 GICs[l-4]. Fig.3a shows Pos(G) of intercalated 2 
and 3L~1623 and 1625cm -1 , very close to that of Stage- 
1 GICs[l|-[3- However, when measuring in-situ the SLG 
Raman spectrum, i.e. keeping the sample sealed in a 
quartz tube after cooling at room temperature following 
a 30 minutes doping, Pos(G) reaches~1627cm -1 , indica- 
tive of heavy doping. This is slightly higher than previ- 
ously reported for Stage- 1 GICs[l|-[4l], and much higher 
than Stage-2 GICs[l|-l4[. This is thus evidence of dop- 
ing on both top and bottom surfaces of our SLG. When 
this sample is removed from the tube and exposed to air, 
Pos(G) goes to~1605cm -1 , similar to Fig. 2, due to des- 
orption. Figs. 3a, b show that Pos(G) for BLG is lower 
than in 3L, and I(2D)/I(G) for BLG stronger than 3L. 
This suggests de-intercalation in BLG flakes. 

To further study the de-adsorption/de- intercalation, 
Figs.4a-c report the BLG Raman spectra as a func- 
tion of time for a period of up to 8 months. Pos(G) 
starts at~ 1623cm -1 , and can be fitted with a single 
Lorentzian with FWHM(G)~8cm -1 , indicating uniform 
doping HI H E3, EH. After twelve hours, Pos(G) 
decreases to~1621cm -1 , and FWHM(G) increases 
to^lOcm -1 . After one month, Pos(G)~1620cm -1 and 
FWHM(G)~12cm -1 . The spectrum evolution indicates 
that FeCl3-intercalated flakes are relatively stable in air 
at room temperature for up to one month. The inter- 
calant Raman modes change little within one month. 
However, after four months, they are not seen anymore, 
while G and 2D acquire a multiple peak profile. This 
confirms that a significant reduction in doping occurs af- 
ter one month due to de-adsorption and de-intercalation. 
The multiple G peaks after 4 months may result from 
different de- adsorption rates on different surfaces. 

Figure 5 plots I(2D)/I(G) and A(2D)/A(G) as a func- 



tion of Pos(G). With increasing Pos(G), i.e. increas- 
ing doping, I(2D)/I(G) and A(2D)/A(G) both decrease. 
I(2D)/I(G) and A(2D)/A(G) of heavily doped SLG are 
close to those of almost fully-doped 2 and 3L, further 
confirming double-surface doping. 

We now consider the dependence of I(2D)/I(G) and 
A(2G)/A(G) on doping and excitation wavelength. We 
use a FeCl3-intercalated Stage- 1 GIC as an example to 
show how to probe the Fermi level by multi-wavelength 
Raman spectroscopy, since in this case the Fermi level is 
well know by independent characterizations [lj. Fig. 6a 
plots the Raman spectra measured at 488, 514, 561, 
593 and 633nm. These are similar to those in Fig. 3. 
For all lasers Pos(G)-1626cm -1 ,FWHM(G)-7cm -1 , in 
good agreement with that previously reported for FeCl3- 
intercalated Stage- 1 GICs[H-[3. At 633nm, the 2D 
peak is almost unobservable, similarly to what reported 
in Refl25l for SLGs doped by bromine. However, in- 
creasing the excitation energy from 2.09eV (593nm) to 
2.54eV (488nm), the 2D peak appears. A(2D)/A(G) and 
I(2D)/I(G) are plotted as a function of excitation energy 
in Fig. 6b. The trend of these intensity ratios can be un- 
derstood considering the full resonant Raman scattering 
process for the 2D band 46|, 48j. Fig. 6c plots the doped 
SLG band structure. For a given laser energy, to acti- 
vate the 2D peak an electron-hole pair must be excited 
in process: a— >>b, and recombined in process c— ^d. These 
transitions differ by the 2D peak energy: 



Et = Ep — huj2D 



(2) 



There are three cases 1) When El and Et are both 
larger than 2Ep, the 2D band can be always observed; 
2) When El is larger than 2 Ep, but Et is smaller than 
2Ep, process c^d is forbidden due to Pauli blocking; 3) 
when both El and Et are smaller than 2Ep, both pro- 
cesses a^b and c— ?>d are forbidden. Therefore, only when 
Et>2Ef, i.e., (El-Tvj02d) /2>Ep 1 the 2D band is observ- 
able. Thus, the absence of the 2D band in the Raman 
spectra of FeCl3-intercalated Stage- 1 GICs, and FeCl3 
and Br 2 heavily-doped SLGsfHI at 1.96eV (632. 8nm) in- 
dicates that their Ep is larger than 0.81eV. When (El- 
hu 2D )/2>E F , both I(2D)/I(G) and A(2D)/A(G) should 
increase. Thus, the sharp intensity increase when mov- 
ing from 2.21eV (561nm) to 2.09eV (593nm), implies that 
Et corresponding to E£=2.09eV is close to 2Ep. Then, 
Ep ~0.88eV, close to that measured in Stage-1 GICs by 
electron energy loss spectroscopy [lj. 

Ep can be also estimated from A(G)/A(2D) at 514nm, 
by using Eq.l. The numerical values of f(e 2 /svp) can 
be taken from Fig. 3 in ReflU and 7e_ p / l =21meV can 
be extracted from the hole-doping data of Ref.13 mea- 
sured at 514nm. While e is not available for intercalated 
FeCl3, we can estimate it from that measured for FeCl3 
in aqueous solutions [51]. This gives /=0.0908. Then, in- 
serting A(2D)/A(G)=0.688 in Eq.l, gives E F -0.84eV, 
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FIG. 6. a) Raman spectra of stage- 1 GIC measured at 488, 
514 532, 561, 593 and 633nm, normalized to have the same 
G band intensity, b) A(2D)/A(G) and I(2D)/I(G) as a func- 
tion of excitation energy, c) Schematic diagram of SLG band 
structure and fully resonant Raman processes for the 2D 
band: a— ?>b, photon absorption; c— )>d, electron-hole recom- 
bination; b— ^c and d— »a, phonon emission. 

very close to that derived by the intensity transition as 
a function of excitation energy. 

In conclusion, graphite flakes consisting of a few 
graphene layers can be doped by adsorption and interca- 
lation of FeCl3. This results in each of the layers behaving 
as a hole-doped SLG. These are stable up to one month 
after air exposure. The variation of the 2D intensity rel- 
ative to the G peak with excitation energy allows one 
to estimate the Fermi energy. We get i?p~0.9eV, corre- 
sponding to a fractional charge transfer of~l/6.6=0.152 
holes per carbon, i.e.~5.8x 10 14 cm -2 , larger than the 
4x 10 14 cm -2 recently reported by employing a solid poly- 
mer electrolyte gate [17]. 
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